Dynamical properties of a nonequilibrium quantum dot close to localized-delocalized 

quantum phase transitions 
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We calculate the dynamical decoherence rate and susceptibility of a nonequilibrium quantum dot 
close to the delocalized-to-localized quantum phase transitions. The setup concerns a resonance- 
level coupled to two spinless fermionic baths with a finite bias voltage and an Ohmic bosonic bath 
representing the dissipative environment. The system is equivalent to an anisotropic Kondo model. 
As the dissipation strength increases, the system at zero temperature and zero bias show quantum 
phase transition between a conducting delocalized phase to an insulating localized phase. Within 
the nonequilibrium functional Renormalization Group (FRG) approach, we address the finite bias 
crossover in dynamical decoherence rate and charge susceptibility close to the phase transition. We 
find the dynamical decoherence rate increases with increasing frequency. In the delocalized phase, 
it shows a singularity at frequencies equal to positive or negative bias voltage. As the system 
crossovers to the localized phase, the decoherence rate at low frequencies get progressively smaller 
and this sharp feature is gradually smeared out, leading to a single linear frequency dependence. 
The dynamical charge susceptibility shows a dip-to-peak crossover across the delocalized-to-localized 
transition. Relevance of our results to the experiments is discussed. 

PACS numbers: 72.15.Qm,7.23.-b,03.65.Yz 



Introduction 

Quantum phase transitions (QPTs)[l|, due to com- 
peting quantum ground states are of fundamental impor- 
tance in condensed matter physics and have attracted 
much attention both theoretically and experimentally. 
Near the transitions, exotic quantum critical properties 
are realized. In recent years, there has been a growing in- 
terest in QPTs in nanosystemsjsHlj]. Very recently, QPTs 
have been extended to nonequilibrium nanosystems with 
a large bias voltage being applied to the setups. Close to 
QPTs much of the attention has been focused on equi- 
librium properties; while relatively less is known on the 
nonequilibrium properties. The key difference between 
equilibrium and nonequilibrium properties near QPTs 
is the voltage-induced nonequilibrium decoherence rate 
which behaves very differently from that in equilibrium 
at finite temperatures, leading to distinct nonequilibrium 
properties near QPTs. 

Recently, two generic examples ll| have been stud- 
ied: (i). the transport through a dissipative resonance- 
level (spinless quantum dot) at a finite bias voltage 
where dissipative bosonic bath (noise) coming from the 
environment in the leads (ii). a spinful quantum 
dot coupled to two interacting Luttinger liquid leads [l2| 
where the electron interactions can be regarded as an 
effective Ohmic dissipative bosonic bath 25]. As dissi- 
pation (or interaction) strength is increased, both sys- 
tems can be mapped onto different effective Kondo mod- 
els, which exhibit QPT in transport from a conducting 



(delocalized) phase where resonant tunneling dominates 
and an insulating (localized) phase where the dissipa- 
tion (or electron-electron interaction) prevails. Similar 
dissipation driven QPTs have been investigated in vari- 
ous systems 0, To obtain the nonequilibrium trans- 
port properties, the authors in Ref. 11 1 and Ref. [l2[ ap- 
plied the nonequilibrium Renormalization Group (RG) 
approach [l^ in the form of self-consistent scaling equa- 
tions for frequency-dependent Kondo couplings and the 
static decoherence rate T(V,T, B). Though the dynam- 
ical nonequilibrium effects in Kondo models have been 
addressed|17i fl8l [20| . less is known of the steady-state 
nonequilibrium decoherence effect on the anisotropic 
and/or two-channel Kondo models. In this paper, we 
address the nonequilibrium decoherence effect on the 
anisotropic Kondo model in the presence of a large bias 
voltage, which is relevant for the delocalized-to-localized 
nonequilibrium QPTs in the dissipative quantum dot 
(the first example mentioned above). To obtain the dy- 
namics (frequency dependence) of decoherence rate, we 
generalize the approach taken in Ref. [Ill Ell Yia a Func- 
tional Renormalization Group (FRG) method developed 
in Ref. [16j |. The nonequilibrium decoherence rate is di- 
rectly proportional to the width of the peak in dynamical 
spin susceptibility. We furthermore investigate the spec- 
tral properties of the dynamical decoherence rate close 
to the QPT and its implications to the dynamical charge 
susceptibility, which can be measured experimentally. In 
particular, as the system goes from the delocalized to the 
localized phase we find the dynamical decoherence rate 




FIG. 1: (Color online) Diagram for the pseudofermion self 
energy. 



for small frequencies gets smaller in magnitude and the 
singular "kink-like" behavior occurring at the frequencies 
equal to the bias voltage (lo — ±V) gets smeared out. We 
have calculated the dynamical spin susceptibility. As the 
system moves from the delocalized to the localized phase, 
it shows a dip-to-peak crossover and the smearing of the 
sharp feature at u) rs ± V. The relevance of our results to 
the experiments is discussed. 



Dissipative resonant level model 

Model Hamiltonian. The system we study here is a 
spin-polarized quantum dot coupled to two Fermi-liquid 
leads subjected to noisy Ohmic environment, which cou- 
pled capacitively to the quantum dot[ll|. For a dissi- 
pative resonant level (spinless quantum dot) model, the 
quantum phase transition separating the conducting and 
insulating phase for the level is solely driven by dissipa- 
tion. The Hamiltonian is given by 11| : 



H = 2J ( e (X) _ Mi) c L c « +ticl i d + h.c. 

k,i=l,2 

+ ^2x r (dU-l/2)(b r + bt) + J2^rbib r , 

r r 

+ h(dU- 1/2) (1) 

where ti is the hopping amplitude between the lead i and 
the quantum dot, Cki and d are electron operators for the 
Fermi-liquid leads and the quantum dot, respectively, 
jttj = ±V/2 is the chemical potential (bias voltage) 
applied on the lead i, while h is the energy level of 
the dot. We assume that the electron spins have been 
polarized by a strong magnetic held. Here, b a are the 
boson operators of the dissipative bath with an ohmic 
spectral density J (uj) = ^2 r \^8(uj ~ uj t ) = auj with 
a being the strength of the dissipative boson bath. 



Through similar bosonization and refermionization 
procedures as in equilibrium the above model is 

maped onto an equivalent anisotropic Kondo model in 
an effective magnetic field h with the effective left L 



and right R Fermi- liquid leads [ll|. The effective Kondo 
model takes the form: 



Hk 



E 

fe,7=i,fl,o-=t4 



fa ~ ^JcLcfc^ 



where c^ L , fl , is the electron operator of the effective 
lead L{R) 1 with spin a. Here, the spin operators 
are related to the electron operators on the dot by: 
S+ = d\ S- = d, and S z = cfid - 1/2 = n d - 1/2 where 
rid = d^d describes the charge occupancy of the level. 
The spin operators for electrons in the effective leads 
are s± p = J2 a ,s,k,k< 1 / 2c l 7 a <T aS c k' 0S, the transverse and 

1(2) 

longitudinal Kondo couplings are given by J ± y ' oc t 1{2) 
and J z oc 1/2(1 — 1/V2a*) respectively, and the effective 
bias voltage is /x 7 = ±y ^/l/(2a*), where 1/a* = 1 + a. 
Note that /i 7 — > ±V/2 near the transition (a* — !• 1/2 
ora-} 1) where the above mapping is exact. The spin 
operator of the quantum dot in the effective Kondo 
model S can also be expressed in terms of spinful 
pseudofermion operator f a : S i=XyVyZ = fX^iix^^fp- 
In the Kondo limit where only the singly occupied 
fcrmion states are physically relevant, a projection 
onto the singly occupied states is necessary in the 
pseudofermion representation [ill 14 1. This can be 



achieved by introducing the Lag range multiplier A so 

In equilibrium, the 



1 13, 15 



that Q = E 7 4/7 
above anisotropic Kondo model exhibits the Kosterlitz- 
Thouless (KT) transition from a delocalized phase with 
a finite conductance G ~ -^-g (e = 7j = 1) for J± + J z > 
to a localized phase for J± + J z < with vanishing 
conductance. The distinct profile in nonequilibrium 
transport near the localized-delocalized KT transition 
has been addressed in Ref. 11]. Below we will turn our 



attention to the dynamical charge susceptibility of the 
quantum dot close to the transition. 

Nonequilibrium FRG formalism. The non-equilibrium 
perturbative FRG approach is based on the generaliza- 
tion of the perturbative RG approach studied in Ref. 1141 
for the nonequilibrium Kondo model. Following Ref. [16j . 
the frequency dependent RG scaling equations for the ef- 
fective Kondo couplings in the Keldysh formulation are 
given bvflij]: 



dh\D 
<91nL> 



E 



pv_ 

2 



/9=-l,l 



h.c. 



where g ± (w) = N(0)Jl = N(Q)Jl„, g z (u) = N(0)J z 
are dimensionless frequency-dependent Kondo couplings 
with N(0) being density of states per spin of the conduc- 
tion electrons (we assume symmetric hopping t\ = t% = 
t). Here, 9 U = 6(L> - \w + iT(u))\) (with D < D„ be- 
ing the running cutoff) comes from the leading logarith- 
mic corrections for the Kondo vertex function originated 
from the product of the Keldysh component of the lead 
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(6) 



> and xtf^J^ reads: 

<(>),o/3. 



Xc£ & = J ^We-tt)(G QCT ( £ + tt)G /CT (fi)) <(>) , 

(7) 

where G is the Green's function in 2 x 2 Keldysh space, 
and its lesser and greater Green's function are related to 
its retarded, advanced, and Keldysh components by: 
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FIG. 2: Frequency dependent Kondo couplings of a dissipa- 
tive resonant level model at zero temperature for (a). g±(ui) 
and (b). g z (uf) across the localized-delocalized transition for 
different bare Kondo couplings via the FRG approach (solid 
lines) and the approach in Ref. [ill ] (dot-dashed lines). We 
have set V = Q.2Dq where Do = 1 for all the figures. 



electron Green function G„ (to) with the real part of the 
retarded or advanced dressed pseudo fermion propagator 



G< = (G K -G R + G A )/2 
G> = (G K + G H -G A )/2 



(8) 



Specifically, the lesser and greater components of Green's 
function of the conduction electron in the leads and of the 
quantum dot (impurity) are given by (in the absence of 
magnetic field): 



Re(Gf /A )[16|: 



'iR I 
T f<r\ 



-2ni tanh(- 



2T 



-)jv e(A) - M) 



GfM = 
GUe) = 



iA c (e){l - fe-p L/R ) 

27ri5(e)n/ CT (e) 

2mS(e)(n fa {e)-l) 



(9) 



GU") 



(4) 



where (w) = Re(T,^(u))) + ^To-(u}) is the impurity self- 



where A c (e) — 2ttNqQ(Dq — e) is the density of states 
of the leads, rif a {e) = j\f a is the occupation num- 
ber of the pseudofermion which obeys rift + n/j. = 1, 
nfcr(e — > 0) = 1/2 in the delocalized phase and rift je — > 
0) -> 0, n/i(e -> 0) -> 1 in the localized phase [ill fl3j|. 
Here, the pseudofermion occupation number rif a and 
the occupation number on the dot rid are related via 
< rift - n/4. >=< n<j > -1/2. Also, f u -p h/R is the 
Fermi function of the L/i? lead given by f u -^ L/a = 



energy (defined below) with the imaginary part being the 
dynamical decoherence rate T a (uj), and Nq = At 
T = 0, the above correction can be approximated by the 
0-function in the RG equations shown above. Also, the 
dynamical decoherence (dephasing) rate T a {uS) at finite 
bias which serves as a cutoff to the RG flow of 5x^(w)Q, VC 1 + ^"-"v*)/** 2 ')). We have therefore 
and it is determined self-consistently along with the RG 
equations for the Kondo couplings. Note that in general 
r (T (w) depends on the impurity spin <r; however, in the 
absence of the magnetic field as we consider here, we 
have the spin symmetry and hence: T^(u)) = r^(w) = 
T(u;). We can obtain T(lu) by the imaginary part of the 
pseudofermion self energy via second-order renormalized 
perturbation theory (see Fig. [TJ: 



Xcj>(<0 = 27r% a (e)A c (e)(l-n /(7 (e))/ 6 _ Mo 



XcT/W = ^igp a (e)A c (e)n fa (e)(l - f e ^ a ) (10) 



r H = 4- / de 9-L(e + w)5_L(e)[/ e _^(l-/ e+w _ MR )] 
+ 5 z (e + w)3 2 (e)[/ e _ ML (l - / e+w _ M J] 



(11) 



10 / 2^ 



-5a/3(^ + e) 



a,P=L,R 

X^ W W<$ (<> (« + «) 



(5) 



where a =j, J,, CT(T ' is the tensor associated with the prod- 



The FRG approach here is accomplished by self- 
consistently solving the RG scaling equation Eq. [3] 
subject to Eq. [5] and Eq. [11] The solutions at zero 
temperature for q \ (uj) and g a}Z (uj) across the transition 
are shown in Fig. [2[ll]. Note that our FRG approach 
is somewhat different from that in Ref. [l6j |: We do not 
formulate and solve for the RG scaling equation for the 
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FIG. 3: (Color online) T(u) (rescaled to g\ + gl with g± :Z 
being the bare Kondo couplings) versus u at zero temperature 
across the delocalized-localized KT transition. Inset: the log- 
log plot of r(w) versus uj. The bias voltage is fixed at V = 
0.2Dq. Here, Do — 1 for all the figures. 



impurity self energy £cr(w) as shown in Ref. [Iff; instead, 
we calculate the imaginary part of the impurity self- 
energy (or the decoherence rate r(w)) self-consistently 
via second-order renormalized perturbation theory. 
Nevertheless, we have checked in the simple Kondo limit 
with isotropic couplings that the frequency-dependent 
renormalized Kondo coupings gx,z{w) obtained here 
agree very well with that obtained via an equivalent 
FRG approach in Ref. [II]). 

As the system goes from the delocalized to localized 
phase transition, the features in g±{uj) at ui — ±V/2 
undergoes a crossover from symmetric double peaks 
to symmetric double dips, while the symmetric two 
peaks in g z (io — ±V/2) still remain peaks. We find the 
above results based on the more rigorous FRG approach 
are in good agreement with the previous heuristic 
method ll| , which provides us with an independent 
and consistency check on the previous results. Note, 



however that from previous approach in Ref. and 
12j the decoherence rate was taken approximately as 
T(lj = 0); we now generalize this T by including the 
frequency dependence. This generalization improves 
the previous RG formalism and it also provides us with 
more features in the dynamical quantities across the 
transition, such as in dynamical charge susceptibility. 
It is worthwhile mentioning that unlike the equilibrium 
RG at finite temperatures where RG flows are cutoff 
by temperature T, here in nonequilibrium the RG 
flows will be cutoff by the decoherence rate T, a much 
lower energy scale than V, T -C V. This explains 
the dips (peaks) structure in g±i z ){ui) in Fig. [2] In 
contrast, the equilibrium RG will lead to approximately 
frequency independent couplings, (or "flat" functions 
g±(ui) « g±. z {u = 0)). In the absence of field h = 0, 



co+v 

FIG. 4: (Color online) Diagram for the charge susceptibility. 
The dressed pseudofermion propagator (double dashed line) 
is calculated via diagram in Fig. [T] 



g±,(z)( U} ) show dips (peaks) at u = ±V/2. We use the 
solutions of the frequency-dependent Kondo couplings 
ff_L,zo-(w) to compute the dynamical charge susceptibility 
of the resonance-level near the transition. 

Dynamical decoherence rate and charge susceptibility. 
We have solved for the dynamical decoherence rate 
r(cj) at zero temperature self-consistently along with 
the RG equations Eq. [31 and the results are shown in 
Fig. [3 As the general trend, we find T(u>) increases with 
increasing frequency; while it decreases in magnitude 
at low frequencies \u>\ < V as the system crossovers 
from the delocalized to the localized phase, is due to In 
additions, in the delocalized phase, it shows a singular 
"kink-like" behavior at frequencies u = ±V, separating 
two different behaviors for \u>\ > V and \ui\ < V. The 
curves of T(ui) for \to\ < V are closer to the linear 
behavior than those for \ui\ > V. As the system moves to 
the localized phase, T(uj) for both \oj\ < V and \u\ > V 
gradually changes its slopes or curvatures and finally 
merges into a single linear behavior deep in the localized 
phase (see Fig. [3j) . 

To understand the above qualitative features, it proves 
useful to simplify the zero temperature dynamical deco- 
herence rate in Eq. [TT] to the following form: 



47T J-V/2-LU/2 J- 



fV/2-u>/2 



}de 



}dc 



-V/2-u)/2 J-V/2+u>/2 

x g ± (e-Lj/2)g ± (e + u/2) 

g r V/2+u>/2 f-V/2+u>/2 
47r JV/2-u/2 J-V/2-U/2 

x g z (e-iu/2)g z (e + uj/2) 1 (\u J \<V), 

r V/2+u>/2 r -V/2+u)/2 

}de 

-V/2-u>/2 JV/2-u>/2 

x g ± (e-Lj/2)g ± (e + uj/2) 

o r V/2+u/2 /■- V/2+w/2 

+ j-i I + w 

47r JV/2-U/2 J-V/2-U/2 

x g 9 {e-w/2)g,(e + w/2),{\u\>V). (12) 



o r V/2+u/2 r - 

FM = j- 1 + 

47r J-V/2-U/2 JV, 



In the "flat" (or in the "equilibrium form") ap- 
proximation where g± tX (u)) are treated as flat func- 
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FIG. 5: (Color online) Xzz( w ) versus u at zero temperature 
across the localized-delocalized KT transition. The bias volt- 
age is fixed at V — 0.2Dq. Here, Do = 1 for all the figures. 
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FIG. 6: (Color online) h(uj) versus uj at zero temperature 
across the localized-delocalized KT transition. The bias volt- 
age is fixed at V = 0.2Do- Here, Do = 1 for all the figures. 



tions g±, z {V 7 uj = 0) ~ g±, z {V = T) with T being 
temperature[ll|, Eq.Q2]can be simplified as the following 
linear behaviors: 

r» = |^i(o) + u^(o)MM<n 
r» = A w ( g i(o) + ^(o)),(M>v) 

(13) 

As shown in Fig. [3j the above approximated form T e (ui) 
agrees well with T(uS) for \uj\ < V, but it shows devi- 
ation from r(w) for \uj\ > V. The approximated form 
T e (uj) exhibits two linear behaviors with different slopes 
for \uj\ < V and > V, respectively, leading to a "kink- 
like" singularity at lo = ±V. As the system moves from 
delocalized to localized phase, the ratio g±(0)/g z (0) be- 
comes progressively smaller, leading to the suppressions 
of the two slopes. Finally, as the system is deeply in the 
localized phase, these two lines merge into a single line 
since |g±(0)| <C |g z (0)| there. The qualitative behaviors 
of r e (u;) can explain the overall monotonically increasing 
trend of T(lj) with increasing the magnitude of frequency 
as well as the decreasing trend of T(\uj\ < V) as the sys- 
tem moves from the delocalized to the localized phase. 
When the full frequency dependence of g±, z (co) is consid- 
ered, we find T{lu) deviates from the perfect linear behav- 
ior in r e (w) (see Fig. [3]). Furthermore, in the delocalized 
phase, the correction to the linear behavior, T(lu) — T e (uj) 
is more noticeable in the large frequency regime (|u;| > V) 
than in the small frequency regime < V (see Fig. [3]). 
This comes as a result of the wider range in energy e 
to be integrated over in Eq. [12] for \u>\ > V, which ac- 
cumulates more deviations from the flat approximation 
due to the dip-peak frequency dependence of gx tX (cj). As 
the system moves to the localized phase, T(\v\ > V) be- 
comes closer to the linear behavior. This is due to the 



fact that g z (oj) becomes flatter in the localized phase, 
giving rise to a more linear behavior for T{oj). Mean- 
while, based on the solutions for g± tZ (io), we find the 
correction to the linear behavior as well as the singulari- 
ties at lo — ±V are logarithmic in nature in the isotropic 
Kondo limit (g± = 0.05 = g z ) and at the KT transition 
(gj_ = 0.05 = — g z ); while they are power-law in nature 
for |<7_l| 7^ \g z \ as g± iZ (ui) show logarithmic and power- 
law properties in these limits, respectively [lH. Note also 
that for \ui\ > V we find T(lj) < r e (w) in the delocalized 
phase; while the opposite holds in the localized phase. 
This can be understood as in the delocalized phase we 
have |<7j_( z )(w)| < |<?j_( z )(0)| for the majority of the fre- 
quencies (except for u> very close to ±V/2), leading to a 
smaller value of T(lu) compared to r e (w); while the op- 
posite is true in the localized phase. We have checked nu- 
merically that T(lu = 0) obtained here indeed reproduces 
the fre quen cy independent decoherence rate T obtained 
inRef. [lj. 

The effect of the dynamical decoherence rate can 
be measured experimentally via dynamical charge sus- 
ceptibility Xc(w): Im(xc(w)) oc lirn/j-^o drid(u))/dh with 
h oc (N — 1/2) being the effective magnetic field mea- 
suring the deviations of N electrons on the dot from 
the charge degeneracy point. Here, Xc(w) is the Fourier- 
transformed charge susceptibility defined as: 

Xc (t) = i6(t) < [(n d (t) - 1/2), M0) - 1/2)] > (14) 

Experimentally, the dynamical charge susceptibility can 
be measured by the capacitance lineshape in an AC field 
near the charge degeneracy point via the high sensi- 
tivity charge sensor in the Single Electron Transistor 
(SET) connected to the dot 2^. The decoherence rate 
T(lo = 0) here corresponds to the broadening of the res- 
onance peak in the imaginary part of Xc(u)- We have 
calculated XcO*0 at zero temperature based on the renor- 



malized second-order perturbation theory (see the dia- 
gram in Fig. [4j. From the mapping mentioned above, 
the dynamical charge susceptibility x c ( w ) is related to 
the z— component of the effective "spin-spin" correlation 
function Xzz m the effective Kondo model through: 

Xc (t) = -iO{t) < {(n d (t) - 1/2), (n d (0) - 1/2)] > 

Xzz (t) = -iO(t) < [S z (t),S z (0)} > (15) 

Taking the Fourier transform of Xzz (t) and evaluating the 
diagram in Fig.21 the imaginary part of X*z(w)> Xzzfa) — 
Ini(Xzz(u)) is given by: 

= £ [G<> + £ )G> ff (e) 

- G> a (u + e)G< a (e)}, (16) 

where G^^ denote the greater (lesser) components of 
the dressed pseudofermion Green's functions: 

Gf a (uj) = 2ninf <r (w)Af ^(w) 
G>» = 2ni(n fa {u) - l)A fir (u) 

A f(T (cj) = Im[ 1 ], (17) 

oj + hf a (LU) + vq 

where rif a (u:) is the nonequilibrium occupation number 
determined self-consistently by the quantum Boltzmann 
equation [lj|: 

n CT H = (l-S>H/E<H)- 1 (18) 

where (w) are defined in Ref. 13]. The resulting 

expression for \zz{oj) at T = is given by: 



r( w ) 



w 2 + r 2 (w) 

where nf^ipS) is given bvflij]: 



x (1 - 2n /t (w)) 



(19) 



n/ t (w) 
n/ t (w) 



ffi(0)(V- W ) 



1 

2 2 s 2(oju + S i(0)(V + w 
l g j(0)(y-c-2 g 2 u;) 
2 ffi(0)(V + w) 



■,(0 < w < V) 



(-V <LO<0). 



(20) 



Note that we have neglected the vertex correction in 
the calculation for % zz (u;) as it gives only a sub- leading 
correction to Eq. [121 

As shown in Fig. [3 in the delocalized phase, Xzz{u — > 
0) oc to — > 0; while Xzz{u — > 0) — > 1/w in the localized 
phase. Hence, as the system crossovers from delocalized 
to localized phase x" zz {w) shows a dip-to-peak crossover 
for small uj. This behavior can be understood as follows. 
In the delocalized (Kondo screened) phase, the effective 
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FIG. 7: (Color online) Xzz( u )/h(u}) (normalized to 1) versus lo 
at zero temperature across the localized-delocalized KT tran- 
sition. Here, h(uj) is defined in the text. Inset: the log-log 
plot of x zz (uj) (normalized to 1) versus cj. The bias voltage 
is fixed at V = 0.2Dq. Here, Do = 1 for all the figures. 



local "spin" get Kondo screened in the low energy scale; 
therefore, the spin susceptibility should vanish. On the 
other hand, in the localized (ferromagnetic) phase, the 
unscreened free "spin" gives rise to the Curie-law sus- 
ceptibility Xzzi^) ^ l/ u at low energies. Meanwhile, 
in the delocalized phase, we find a a "kink-like" singu- 
lar behavior in x zz (w) at u) = V, coming from the sin- 
gular behaviors of both T(td = V) (see Eq. [3]) and the 
factor 1 — 2n/f in x zz ( w — ^0 ( see Eq. US and discus- 
sions below). However, this singularity gets smeared out 
as the system crossovers to the localized phase. We have 
checked that our results for x zz (^) in the isotropic Kondo 
limit agrees qualitatively well with those in Ref. and 
Ref. [19(. For u ^> V, we find Curie-like susceptibility 
X zz (w) cx 1/lu in both localized and delocalized phases, 
followed Eq. [TJ and Eq. [T3 

We can furthermore link our results to the equilibrium 
and nonequilibrium fluctuation-dissipation theorem (l7l- 
21] . It is useful to define the dynamical fluctuation- 
dissipation ratio h(uj): 



Szzi^y 



(21) 



where S zz {uj) is the Fourier-transformed longitudinal 
spin-spin correlation function with its real-time form 
given by: 



S zz (t) = ±<{S z (t),S z (0)}> 



(22) 



The dynamical spin-spin correlation function S(uj) is 
given by: 



5 M = ft £ ^> + e ) 6 />) 

a=T4 
=<> I, , \ r \r<< 



+ G>> + e)G<»], 



(23) 
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Carrying out similar calculations as in x^(w), we find 
the fluctuation-dissipation ratio h{oj) reads: 

h(u) = 1 - 2n /t (w), (24) 

where n^{ui) is defined in Eq. 1201 In equilibrium, the 
ratio h(uj) = h°(oj) respects the fluctuation-dissipation 
theorem [l7H2C| . given by h°(u>) = 1 — 2n*(u;) = 
tanh(/3w/2) where n®{uj) = eA ii +1 is the Fermi function 
for pseudofermion in equilibrium. At T = (/3 — > oo), 
we have h (u> > 0) = 1, the signature of the equilibrium 
fluctuation-dissipation theorem at T = 0. In nonequi- 
librium and at T = 0, however, we find in general a 
deviation of h(w) in the delocalized phase from the equi- 
librium fluctuation-dissipation theorem: h(u>) < 1 for 
u < F[H Sq|, and h(u -> 0) cx w -> (see Fig. [B] 
and Eq. [Ml Eq. |2"0)) : while we recover the equilibrium 
fluctuation-dissipation theorem h(u) = 1 for M > V, in 
agreement with the results in found in Ref. [171 ]. How- 
ever, that as the system moves to the deep localized 
phase, the above deviation gets smaller and finally we 
recover the equilibrium fluctuation-dissipation theorem 
for all frequencies: h{uS) = 1 for all \cj\ < Do (see 
Fig. [6]) . To more clearly observe the singular property of 
X (M = V), we also plot the normalized Lorentzian form 
r(Q)x"(w)/fc(w) = T(u)T(0)/(lo 2 +r 2 (w)) (see Fig. El). 
As the system moves from the delocalized to the localized 
phase, the width of the Lorentzian peak gets narrower 
and the singularity at cu = V gradually disappears. 

Conclusions 

In conclusion, we have calculated the zero temperature 
nonequilibrium dynamical decoherence rate and charge 
susceptibility of a dissipative quantum dot. The sys- 
tem corresponds to a nonequilibrium anisotropic Kondo 
model. We generalized previous RG approach based on 
Anderson's poor man's scaling approach and renormal- 
ized perturbation theory to a more rigorous functional 
RG (FRG) approach. Within the FRG approach, as the 
systems crossover from the delocalized to the localized 
phase, we find a crossover in dynamical charge suscep- 
tibility from a dip to a peak structure for \u\ < V. 
Meanwhile, we find a smeared-out of the singular be- 
havior in charge susceptibility at u> = ±V in the above 
crossover. We also show the deviation of the equilib- 
rium fluctuation-dissipation theorem for small frequen- 
cies u> < V; while the theorem is respected when the sys- 
tem is in the extremely localized phase. The above sig- 
natures can be used to identify experimentally the above 
localized-delocalized quantum phase transition out of 
equilibrium in nano-systems associated with the Kondo 
effect. 
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